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We present a variational calculation of the spin wave excitation spectrum of double-exchange 
ferromagnets in different dimensions. Our theory recovers the Random Phase approximation and 
f/S expansion results as limiting cases and can be used to study the intermediate exchange coupling 
and electron concentration regime relevant to the manganites. In particular, we treat exactly the 
long range three-body correlations between a Fermi sea electron-hole pair and a magnon excitation 
and show that they strongly affect the spin dynamics in the parameter range relevant to experiments 
in the manganites. The manifestations of these correlations are many-fold. We demonstrate that 
they significantly change the ferromagnetic phase boundary. In addition to a decrease in the magnon 
stiffness, we obtain an instability of the ferromagnetic state against spin wave excitations close to the 
Brillouin zone boundary. Within a range of intermediate concentrations, we find a strong softening 
of the spin wave dispersion as compared to the Heisenberg ferromagnet with the same stiffness, 
which changes into hardening for other concentrations. We discuss the relevance of these results to 
experiments in colossal magnetoresistance ferromagnets. 
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I. INTRODUCTION AND PROBLEM SETUP 

The interactions between itinerant carriers and local 
magnetic moments lead to new magnetic properties in 
a wide variety of systems that have been the subject of 
intense research lately. Examples where such magnetic 
exchange and Hund's rule interactions play an impo- 
rant role range from ferromagnetic semiconductors such 
as EuO, EuS, chrome spinels, or pyrochlorei to dilute 
III-Mn-V and II- VI magnetic semiconductors!^ Of par- 
ticular interest here are the manganese oxides (man- 
ganites) Ri-zA-rMnOa, where R=La, Pr, Nd, Sm, • • • 
and A= Ca, Ba, Sr, Pb, • • • , which display colossal 
magnetoresistancei&i The above systems display ferro- 
magnetic order mediated by itinerant carriers. Our main 
goal in this paper is to describe the role of ubiquitous 
three-body correlations (beyond the mean field approxi- 
mation) on the spin dynamics of such ferromagnets. 

Given the wide variety of such ferromagnetic systems, 
it is important to understand the properties of the min- 
imal Hamiltonian that describes the common properties 
and spin dynamics. In order to calculate the effects of 
correlations beyond the mean field approximation, it is 
often necessary to neglect particularities of the individ- 
ual systems, such as chemical structure and crystal en- 
vironment. The most basic model that applies to all 
such materials is the Kondo lattice or double exchange 
Hamiltonian H = K + H exc h + H super + Hjj, where 

K = J2ka SkC^Cko- is the kinetic energy of the itinerant 
carriers. In the manganites (Ri-zA^MnOs), n = 1 — x 
itinerant electrons per Mn atom occupy the band of Mn 
d-states with e g symmetry. We simplify the calcula- 
tion of the three-body correlations of interest here by 
considering a single tight-binding band of cubic symme- 
try and neglect the bandstructure and the degeneracy 



of the e g states (which may be lifted by external per- 
turbations). The operator cj^ creates an electron with 

momentum k, spin ex, and energy = — 2t J3*=i cos fcja, 
where d = 1,2,3 is the system dimensionality and a is 
the lattice constant. From now on we take a = 1 and 
measure the momenta in units of h/a. 

A common feature of all the systems of interest here 
is the strong magnetic exchange interaction, H exc h, be- 
tween the itinerant carrier spins and the local spin-S 
magnetic moments Si, which are located at the N d lat- 
tice sites Hi. In the manganites, these S=3/2 spins are 
due to the three electrons in the tightly bound t2 S or- 
bitals. Introducing the collective localized spin opera- 
tor S q = 1/y/N Sje~ l<lRj and the corresponding spin 
lowering operator S~ = S* — iS^, we express the local 
magnetic exchange interaction in momentum space: 

V kqcr 
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where a = ±1. In the manganites, J > describes 
the ferromagnetic Hund's rule coupling between the lo- 
cal and itinerant spins on each lattice site. H super is the 
weak antiferromagnetic direct superexchange interaction 
between the spins localized in neighboring sites, while H\j 
describes the local Coulomb (Hubbard) repulsion among 
the itinerant electrons. The precise values of the pa- 
rameters entering in the above Hamiltonian are hard to 
calculate for strongly coupled many-body systems such 
as the manganites. Although the parameter estimates 
vary in the literature, typical values are t ~ 0.2-0.5 eV 
and J ~ 2eV, which corresponds to 4 < J/t < 10 i On 
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the other hand, the antiferromagnetic superexchange in- 
teraction is weak, ~ O.Olt. The electron concentration, 
n = (N e /N) d = 1—x where N e is the number of electrons, 
varies from to 1. Ferromagnetism in the metallic state 
is observed within a concentration range 0.5 < n < 0.8 in 
both 3D and quasi-2D (layered) systems. In this paper 
we neglect for simplicity the effects of H super and Hjj (to 
be studied elsewhere) in order to focus H exc h- 

Given the large values of J/t in most systems of in- 
terest, a widely used approximation is the J — ► oo limit 
(double exchange ferromagnet) & In this strong coupling 
limit, the itinerant carrier is allowed to hop on a site only 
if its spin is parallel to the local spin on that site. The ki- 
netic energy is then reduced when all itinerant and local 
spins are parallel, which favors a ferromagnetic ground 
state (double exchange mechanism) . We denote this fully 
polarized half-metallic state by \F) and note that it is 
an eigenstate of our Hamiltonian H . This state describes 
local spins with S z = S on all lattice sites and a Fermi 
sea of spin-| itinerant electrons occupying all momentum 
states with £k < Ef, where Ef is the Fermi energy. 

Another commonly used approximation is to treat 
the local spins as classical (S — > oo limit) A The ferro- 
magnetism can then be described by an effective near- 
est neighbor Heisenberg model with ferromagnetic in- 
teraction. The quantum effects are often taken into 
account perturbatively in 1/S. This 1/S expansion can 
be implemented systematically by using the Holstein- 
Primakoff bosonization methodi2*iS To 0(1/5*), this 
method gives noninteracting Random Phase Approxi- 
mation (RPA) magnons, whose dispersion is determined 
by the exchange interaction>iiii2i In the strong coupling 
limit J/t — > oo, this RPA dispersion coincides with that 
of the nearest neighbor Heisenberg ferromagnet. The 
0(1/ S 2 ) correction to the spin wave dispersion how- 
ever deviates from this Heisenberg form. This correc- 
tion comes from the scattering of the RPA magnon with 
the spin— f electron Fermi sea when treated to lowest or- 
der in the electron-magnon interaction strength (Born 
approximation) AiS 

The role of nonperturbative carrier-magnon correla- 
tions (beyond 0(1/ S 2 )) has been studied by exact di- 
agonalization of small and ID systemsi2ii4 or by us- 
ing variational wavefunctions 1618 inspired from the Hub- 
bard model and the Gutzwiller wavefunctionii2i22i2ii22 
The variational calculations of Refs ll^ll^ treat the lo- 
cal correlations expected to dominate in the strong cou- 
pling limit. 2 — The ferromagnetic (Nagaoka) state \F) was 
shown to become unstable with increasing electron con- 
centration due to the softening of either single particle 
spin excitations^i or long wavelength spin wave excita- 
tions (negative stiffness) The spin wave dispersion 
deviated from the Heisenberg form for very small elec- 
tron concentrations^ In the concentration range 0.5 < 
n < 0.8 relevant to the manganites, Wurth et.alw& found 
small deviations from the Heisenberg form. A similar 
conclusion was reached based on the 1/S expansion^ For 
n ~ 0.7, the 0(1/ S 2 ) magnon dispersion showed a rela- 



tive hardening at the zone boundary in the strong cou- 
pling limit. 10 Given the questions raised in the literature 
about the adequacy of the simple double exchange model 
for explaining the magnetic and transport properties of 
the manganites 7 , it is important to treat the Hamiltonian 
H in a controlled way in the parameter regime relevant to 
the manganites. Such a treatment would allow us to as- 
sess the accuracy of the commonly used approximations 
and understand the successes and limitations of the very 
basic model in explaining the experiments. 

Due to the interplay between the spin and charge de- 
grees of freedom, a good understanding of the spin dy- 
namics is important for understanding the physics of 
colossal magnetoresistance and transport in the man- 
ganites. Several experimental studies of the spin wave 
excitation spectrum have been reported in the litera- 
ture. Heisenberg-like magnons were observed in the fer- 
romagnetic regime for high electron concentrations (typ- 
ically n > 0.7). 23 This observation is consistent with 
the spin wave spectrum obtained in the classical spin 
limit or by using the strong coupling RPA. However, 
for lower electron concentrations 0.5 < n < 0.7, unex- 
pectedly strong deviations from the short range Heisen- 
berg magnon dispersion were observed in several differ- 
ent manganitesi 24 i 25 i 26 i 27 i 28 i 29 i 30 Most striking is the pro- 
nounced softening of the spin wave dispersion and short 
magnon lifetime close to the zone boundary, which indi- 
cate a new spin dynamics in the metallic ferromagnetic 
phase for intermediate electron concentrations 0.5 < n < 
0.7. The physical origin of this dynamics remains un- 
der debate. It has been conjectured that the coupling to 
additional degrees of freedom not included in the dou- 
ble exchange Hamiltonian H is responsible for this new 
spin dynamics. Some of the mechanisms that have been 
proposed involve the orbital degrees of freedom, the spin- 
lattice interaction, the local Hubbard interaction, band- 
structure effects, etci 9 i 28 i 31 i 32 i 33 . 

In this paper we study variationally the low energy spin 
excitations of the half-metallic fully polarized state \F). 
Our focus is on the role of correlations. Our theory treats 
exactly up to three-body correlations between a magnon 
and a Fermi sea pair by using the most general varia- 
tional wavefunction that includes up to one Fermi sea 
pair excitations. As already noted in the context of the 
Hubbard modelfSS* 2 ^ the Gutzwiller wavefunction, which 
treats local correlations ji£ii£ is a special case of such a 
wavefunction. Here we treat both local and long-range 
correlations on equal footing in momentum space in or- 
der to interpolate between the weak and strong coupling 
limits with the same formalism. We treat nonperturba- 
tively in a variational way the multiple electron-magnon 
and hole-magnon scattering processes that lead to vertex 
corrections of the carrier-magnon interaction. The above 
two scattering channels are coupled by three-body corre- 
lations. We show that this coupling is important for the 
intermediate electron concentrations and exchange inter- 
actions relevant to the manganites, while for small (large) 
n the electron-magnon (hole-magnon) scattering chan- 
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nel dominates. In the case of the ID Hubbard model, 
a similar three-body treatment gave excellent agreement 
with the exact results^ Analogous calculations were per- 
formed to describe the electron-Fermi sea pair local Hub- 
bard interactions22i22iM and the valence (or core) hole- 
Fermi sea pair interactions that lead to the Fermi Edge 
( X-ray Edge) Singularityi^^ 

Our variational wavefunction offers several advantages. 
While local correlations^ 6 *^ dominate in the strong cou- 
pling limit, long range correlations become important 
as J/t decreases^ By working in momentum space, we 
treat both long and short range correlations while ad- 
dressing both the weak and strong coupling limits with 
the same formalism. We therefore expect that our re- 
sults interpolate well for the intermediate values of J/t 
relevant to the manganites Our wavefunction satisfies 
momentum conservation automatically, which reduces 
the number of independent variational parameters. Fur- 
thermore, our results become exact in the two limits of 
N e = 1 and N e = N d . We therefore expect that they 
interpolate well for intermediate electron concentrations 
< n < 1. Our variational equations contain the RPA 
(0(1/5*)), ladder diagram approximation, and 0(1/S 2 ) 
results as special cases. Finally, our results converge with 
increasing system size N and thus apply to the thermody- 
namic limit. The only restriction is that we neglect con- 
tributions from two or more Fermi sea pair excitations. 
Such multipair contributions are however suppressed for 
large S, while their contribution in the case of the ID 
Hubbard model was shown to be small. 34 36 

Here we address a number of issues regarding the ef- 
fects of up to three-body carrier-magnon correlations on 
the spin dynamics predicted by the simplest double ex- 
change Hamiltonian. First, by comparing to the 1/S 
expansion, RPA, and ladder approximation results, we 
show that vertex corrections and long range three-body 
magnon-Fermi sea pair correlations, which couple the 
electron-magnon and hole magnon scattering channels, 
play an important role on the spin dynamics in the pa- 
rameter regime relevant to the manganites. We find large 
deviations from the strong coupling double exchange spin 
wave dispersion, including a strong magnon softening at 
the zone boundary in the intermediate electron concen- 
tration and exchange interaction regime. On the other 
hand, for small electron concentrations (relevant e.g. in 
III(Mn)V semiconductors 4 * 5 -) , the electron-magnon mul- 
tiple scattering processes dominate. However, even in 
this regime the deviations from the RPA and 0(1/ S 2 ) 
magnon dispersions can be strong. 

Second, by using an unbiased variational wavefunc- 
tion, we determine the change in the ferromagnetic phase 
boundary due to the three-body correlations and carrier- 
magnon vertex corrections (not included to 0(l/5 2 )). 
The variational nature of our calculation allows us to rig- 
orously conclude that the ferromagnetic state \F) is un- 
stable when its energy exceeds that of the variational spin 
wave energy. In addition to the long wavelength softening 
and eventual instability, which occurs in all dimensions, 



we find another instability for momenta close to the zone 
boundary while the stiffness remains positive. This insta- 
bility only occurs in 2D and 3D for intermediate electron 
concentrations (0.4 < n < 0.7 for the 2D three-body cal- 
culation). This effect is exacerbated by the three-body 
correlations. One should contrast the above instability to 
the spin wave softening (but not instability) at the zone 
boundary that occurs for small n < 0.3M^1& 

Third, we study the deviations from the Heisenberg 
spin wave dispersion induced by the three-body correla- 
tions and vertex corrections. This comparison is impor- 
tant given the experimental observation of pronounced 
deviations for n < 0.7 i 24 i 25 i 26 ' 2 7i 2 9i 2 ?i 3( i ) Deviations from 
Heisenberg behavior already occur to 0(1/ S 2 ), or even 
to 0(1/5*) for finite J/t, but in most cases correspond 
to magnon hardening^ By comparing our results to the 
Heisenberg dispersion with the same stiffness, we show 
that, for values of J/t relevant to the manganites and 
such that the ferromagnetic state is stable up to n ~ 0.8 
or higher, the three-body correlations in the 2D system 
give magnon hardening at the zone boundary for n < 0.4 
followed by strong magnon softening for 0.4 < n < 0.7 
and then small magnon hardening for n > 0.7. This be- 
havior is similar to the experiment. 

The outline of this paper is as follows. In Section II 
we discuss the four approximations that we use to cal- 
culate the effects of the carrier-magnon correlations on 
the spin wave dispersion. In Section II. A we discuss the 
variational wavefunction that treats the three-body cor- 
relations. The variational equations that determine the 
spin wave dispersion are presented in Appendix I. We also 
obtain the RPA magnon dispersion variationally. In Sec- 
tion II. B we establish the connection between the above 
variational results and the 1/S expansion resultsi^iiii We 
show, in particular, that the 0(l/5 2 ) magnon disper- 
sion can be obtained from our variational equations by 
treating the carrier-magnon scattering to lowest order in 
the corresponding interaction strength (Born approxima- 
tion). In Section II. C we discuss the two-body ladder ap- 
proximation, obtained from our variational results by ne- 
glecting the coupling between the electron-magnon and 
hole-magnon scattering channels. The latter coupling 
is discussed further in Appendix II. In Section II. D we 
discuss the approximation of carrier-localized spin scat- 
tering and show that this variational treatment improves 
on the RPA while making the numerical calculation of 
three-body effects feasible in much larger systems. In 
Section III we present our numerical results for the spin 
wave dispersion, ferromagnetic phase diagram, and devi- 
ations from Heisenberg dispersion in the ID, 2D, and 3D 
systems and compare between the different approxima- 
tions. We end with the conclusions. 



II. CALCULATIONS 

In this section we discuss the four approximations that 
we use to treat the effects of the carrier-magnon correla- 
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tions. From now on we measure the energies loq of the 
spin wave states with respect to that of the fully polar- 
ized (half-metallic) ferromagnetic state \F), whose sta- 
bility and low energy spin excitations we wish to study. 
We note that \F) is an exact eigenstate of the Hamil- 
tonian H with maximum spin value and total spin z- 
component N(S + n/2). To study the stability of this 
state in a controlled way, it is important to use approxi- 
mations that allow us to draw definite conclusions. This 
is possible with the variational principle, which allows us 
to conclude that a negative excitation energy ujq means 
instability of \F), driven by the spin wave of momentum 
Q. Our variational states have the form |Q) = Mq\F), 

where the operator Mq conserves the total momentum, 
lowers the z-component of the total spin by 1, and in- 
cludes up to one Fermi sea pair excitations. A spin wave 
has total spin z-component of N(S + n/2) — 1, which cor- 
responds to one reversed spin as compared to \F). This 
spin reversal can be achieved either by lowering the local- 
ized spin z-component by 1 or by coherently promoting 
an electron from the spin— "f band to the spin— J, band. 
The spin reversal can be accompanied by the scattering 
(shakeup) of Fermi sea pairs. From now on we use the 
indices v, ji, ■ ■ ■ to denote single electron states inside the 
Fermi surface and a, (3, • ■ • to denote states outside the 
Fermi surface. 



A. Three— body Correlations 

In this section we discuss our three-body variational 
calculation of the spin wave dispersion. First, however, 
we show that the well known RPA magnon dispersionii*i2i 
can be obtained variationally for any value of J/t by ne- 
glecting in Mq all Fermi sea pair excitations. The most 
general such operator has the form 



M, 



QRPA 



i<V T . (2) 



where the variational equations for the N e amplitudes 
X PPA are obtained in Appendix I (Ea. (|A4| ). In the 
strong coupling limit JS — > oo, X RPA — > 1 and Mq RPA 
reduces to the total spin operator. To lowest order in 
t/JS we obtain from Eqs. l|A3jl and l|A4jl that the RPA 
dispersion then reduces to the Heisenberg dispersion 



,RPA 



1 



2N S 



— Y 

I- n/2 ^ 



(e v „ +Q -£„) + 0{t/JS). (3) 



The 0(1/5*) magnon dispersionifiiil is obtained from the 
above strong coupling RPA result by replacing the total 
spin prefactor S + n/2 by S. 

We now include in Mq the most general contribution 



of the one Fermi sea pair states: 
MX = 



T c mT 



1 



2VN 



CI; fl. 



,(4) 



where the amplitudes X®, \I>Q and $^ are all deter- 
mined variationally; we do not use the RPA results for 
and As compared to previous calculations, we 
do not assume any particular form or momentum depen- 
dence for the above variational amplitudes. This allows 
us to treat in an unbiased way the long range correlations 
for any value of J/t. The first two terms on the rhs of 
Eq.Q create a magnon of momentum Q. The last two 
terms describe the scattering of a momentum Q magnon 
to momentum Q + v — a with the simultaneous scatter- 
ing of a Fermi sea electron from momentum u < Uf to 
momentum a > kp. 

Our wavefunction Eq. (J2J) becomes exact in the two 
limits of N e = 1 and N e = N d . To see this, we note that, 
for N e = 1, the Fermi sea consists of a single electron. 
As a result, multipair excitations do not contribute, while 
$ Q = 0. In the half-filling limit N e = N d , all lattice sites 
are occupied by one spin-"f electron and the Fermi sea oc- 
cupies all momentum states up to the zone boundary. As 
a result, the RPA wavefunction Eq. J2J becomes exact. 
Eq. J3J also gives the exact wavefunction in the atomic 
limit t — 0, Sk — 0, where the variational amplitudes do 
not depend on the electron momenta. To see this, we 
note that, due to the Pauli principle, must be anti- 

symmetric with respect to the exchange of the Fermi sea 
electron momenta v and /x. In the atomic limit, $^ must 
therefore vanish since it is independent of the momenta. 
For the same reason, all multipair amplitudes vanish as 
well and Eq.(0} gives the exact result. 

The variational equation for is derived in Appendix 
I fEas. (|A5|l and l|A8(l ). The magnon energy is obtained 
from Eq. (|A1|) after substituting Eq. I|A2() : 



2N ^JS + £ 



— e v — (j, 



Q 



where we introduced the electron vertex function 



rr 



2N JS- 



- £iH-Q 



wq 



/ j ol'V 



(5) 



(6) 



The first term in Eq. (|5jl gives the RPA contribution to the 
magnon energy. The second term is the carrier-magnon 
self energy contribution, determined by the electron ver- 
tex function T e . The latter satisfies Ea. (|B3p . which de- 
scribes the multiple electron-magnon scattering contri- 
bution (ladder diagrams, two-body correlations) as well 
as the coupling to the hole vertex function, Eq. (|B2Jl . 
due to the three-body correlations. Below we discuss 
three contributions to the full r e : 0(1/ S 2 ) (Born scatter- 
ing approximation), ladder diagram (two-body carrier- 
magnon correlations) , and the contribution due to carrier 
scattering with the localized spins. 
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B. 1/S Expansion 

In this section wc make the connection with the 
Holstein-Primakoff bosonization treatment of the quan- 
tum effectsAiS In particular, we show that the 0(1/ S 2 ) 
magnon dispersion 9 ^ can be obtained by solving 
Eq. (IB3|) perturbatively, to lowest order in the carrier- 
magnon interaction and 1/S. 

First, we recall that classical spin behavior is obtained 
in the limit S — * oo with J5=finite. By expanding Eqs 
(|A8I) . © and l|B2|) in powers of the small parameter 1/S 
(JS= finite) we see that #Q = 0(1/ S), T e = 0(1/S 2 ), 
and r' 1 = 0(1/S 2 ). In particular, we obtain from 
Ea. (jB3|) to lowest order in 1/S 



r e = 



J 



AN 2 \ JS + e v+ Q-e v 



E 



i 



(7) 



and from Eq.JSJ 
J 



E 



2N ^ JS + £„+q - e v - uj, 



J 



E 



4iV 2 tlT (JS + e v+Q - e v f (e a - e v ) 



(8) 



The last term in the above equation comes from the low- 
est order magnon-electron scattering contribution. The 
0(1/ S) spin wave dispersion^ is obtained from the first 
term by neglecting ujq = 0(1/ S) in the denominator. 
The spin wave energy to 0(1/ S 2 ) is obtained by expand- 
ing the first term to this order. We recover the strong 
coupling 0(1/ S 2 ) results of Refs l^fiol obtained by using 
the bosonization technique, by further expanding Eq.© 
in the limit JS — > oo. 

The 0(1/ S 2 ) magnon dispersion is not variational. 
Thus we cannot definetely conclude instability of the fer- 
romagnetic state if we find a negative magnon energy to 
0(1/ S 2 ). On the other hand, the three-body calcula- 
tion outlined in the previous section treats the magnon- 
Fermi sea pair interaction variationally rather than per- 
turbatively (as in the 1/S expansion) while recovering 
the 0(1/ S 2 ) results as a special case. The n-pair con- 
tributions to Eq.© have amplitudes of order 0(1/ S n ). 
Therefore, the shake-up of multipair excitations is sup- 
pressed for large S. Our three-body calculation thus 
puts the 0(1/ S 2 ) results on a more quantitative (varia- 
tional) basis by treating fully rather than perturbatively 
all contributions of the one Fermi sea pair states. 

adaa 



C. Two— body Ladder Approximation 

To go beyond the Born approximation (0(1/ S 2 )), we 
first consider the two-body correlation contributions to 
the Fermi sea pair amplitude while still neglecting the 
three-body correlations. This is equivalent to treating 



the ladder diagrams that describe the multiple electron- 
magnon and hole-magnon scattering, while neglecting 
the coupling between these two scattering channels. Not- 
ing that the magnon dispersion is determined by T e only 
(Eq.©), the ladder approximation dispersion is obtained 
from Eq. {B3J| with T h = and Eq.©: 



J 



J 



AN 2 ^ i _ 



2N ^ JS + e„ + Q - e v -uq 

( £i;+q-g^-"Q \ y i/aQ 



2N JS+e„ + Q-£„-w Q 



2—Y ,1/A Q , 



(9) 



where is given by Ea. ljA7|) . We note that, similar to 
the 1/S expansion, the above ladder approximation re- 
sult is not variational. A similar approximation was used 
in the context of the Fermi Edge Singularityi2£i2& There 
it was shown that at least three-body correlations are 
necessary in order to describe the unbinding of the dis- 
crete exciton bound statei 35 ' 36 In the case of the Hubbard 
model, the ladder approximation was shown to overesti- 
mate the electron self energy^ 

The difference between the spin wave dispersions 
Eq.© and the full three-body calculation comes from 
the three-body correlations. By comparing the corre- 
sponding curves in the next section we can therefore 
judge the role of these correlations on the spin dynamics. 



D. Carrier— localized spin scattering ($ = 0, \P 7^ 0) 

To describe the three-body correlations, the coupled 
equations for T e and T h must be solved. Although this 
is possible in ID and 2D for fairly large systems, in 3D 
the numerical solution of the full variational equations 
is challenging, due to the dependence of on six mo- 
mentum components. On the other hand, depends 
on one momentum only. The dependence of T^ v on the 
momentum v can be eliminated by considering a simpler 
variational wavefunction, obtained from Eq.© by setting 
<f> = 0. This corresponds to treating fully the scattering 
of the electron with the localized spins while neglecting 
the electronic contribution to the scattered magnon. This 
approximation becomes exact in the two limits N e = 1 
and N e = N d and also recovers the 0(1/ S 2 ) and RPA 
results. Its main advantage is that it improves the RPA 
by allowing us to treat variationally three-body carrier- 
localized spin correlations in a large system. The corre- 
sponding spin wave dispersion is obtained by solving the 
coupled Eas. (|B3|) and HB4(1 and then substituting T e in 
Eq.©. 



III. NUMERICAL RESULTS 

In this section we present the results of our numer- 
ical calculations. First of all, we compare our numer- 
ical results with results made by exact diagonalization 




Q Q 

FIG. 1: Spin-wave dispersion at the J — + oo limit in a ID 

system with N =20, S =1/2 and for number of electrons Ne FIG. 2: Spin-wave dispersion in the ID system: comparison 
=1, 3, 5, 7. The points in the figure show exact results. of the full three-body variational calculation (solid curve) to 

the different approximations discussed in the text. J/t = 10. 



of small ID systemic*. Fig 1 shows the comparison 
where we see that our results are in very good agree- 
ment with the exact calculations. To draw conclusions 
on the role of the correlations, we compare the differ- 
ent approximations discussed in the previous section for 
a d-dimensional lattice with N d sites. We perform our 
calculations for d=l, 2, and 3. The dimensionality of 
the system affects the quantum fluctuations and cor- 
relation effects. Quantum fluctuations are expected to 
be most pronounced in the ID system, where we show 
that the 1/S expansion can lead to spurious features. 
The calculation of the ID magnon dispersion could also 
be relevant to quasi-lD materials with chain structures. 
Our 2D magnon dispersion is relevant to the quasi-2D 
layered manganites, 37,3 — where a pronounced spin wave 
softening and deviations from the Heisenberg disper- 
sion similar to the 3D system2i2£i2§i2L2& were observed 
experimentally. 29 - 30 The similarity of the spin dynamics 
in the 3D and 2D systems indicates that the relevant 
physical mechanisms are generic and do not depend cru- 
cially on the particularities of the individual systems. In 
2D, the full three-body variational calculation can be 
performed in fairly large systems ( N ~ 20 — 30) , while in 
3D it could only be performed for N ~ 10. Therefore, the 
2D system also offers computational advantages. On the 
other hand, the rest of the approximations discussed here 
can be performed in very large systems (up to N ~ 200) , 
until full convergence with increasing N is reached. 

We start with the dependence of the spin wave excita- 
tion spectrum on the electron concentration n. Figs. 2 
and 3 show the magnon dispersion in the ID and 2D 
systems respectively for a fixed exchange intreraction, 
J/t = 10, and four different values of n. The 2D dis- 
persion (Fig. 3) was calculated along the Brillouin zone 
direction (0, 0) — > (ir, 0) (T — X), where the discrepancies 
between the different approximations are maximized. 

For very small electron concentrations [n = 0.12 in 



Figs. 2 and 3), the carrier-magnon scattering tends to 
soften the spin wave dispersion close to the zone bound- 
ary, consistent with previous results^S^ This can be 
seen by comparing in Figs. 2(a) and 3(a) the RPA disper- 
sion with the different calculations of the carrier-magnon 
scattering effects. It is important to note that, despite 
the softening, the spin wave dispersion does not become 
negative (unstable) close to the zone boundary for very 
low concentrations. This softening may be interpreted 
as a remnant in the thermodynamic limit of the failure 
of the RPA for N e = 1, where Eq. J3J gives the exact 
solution. Indeed, for N e — 1, the magnon energy is of 
0(1/N 2 ), while the RPA gives 0(1/N) energies. Fur- 
thermore, for low concentrations, the $ 7^ 0, \& ^ and 
<I> = 0, W 7^ variational calculations give results similar 
to the ladder approximation (the corresponding curves 
almost overlap in Figs. 2(a) and 3(a)). This indicates 
that the three-body correlations are weak for very low 
concentrations. This result can be understood by noting 
that the last term in Eq. HA8fl (and Eq. jA70). which 
describes the hole-magnon multiple scattering contribu- 
tion, is suppressed for small n. Indeed, with decreas- 
ing n and Fermi energy Ep, the phase space available 
for the hole to scatter decreases relative to the phase 
space available for electron scattering. As a result, the 
electron-magnon scattering channel (electron ladder di- 
agrams) dominates. On the other hand, the difference 
between the above dispersions and the RPA is strong, 
while the differences from the 0(1/ S 2 ) (Born scattering) 
result are noticable even for very small n (Figs. 2(a) 
and 3(a)). The latter differences come from the multiple 
electron-magnon and hole-magnon scattering processes. 
In ID, the 0(1/ S 2 ) result fails qualitatively for very low 
(n = 0.12 in Fig. 2) and very high (n > 0.8) electron con- 
centrations. For such concentrations, the 0(1/ S 2 ) dis- 
persion becomes negative (unstable) at the zone bondary. 
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FIG. 3: Spin- wave dispersion in the 2D system along the 
direction F — X for the same parameters as in Fig. 2: com- 
parison of the full three-body variational calculation (solid 
curve) to the different approximations discussed in the text. 



This zone boundary instability persists even in the strong 
coupling limit J — * oo but is absent in all our variational 
results. 

With increasing electron concentration, the spin wave 
energies and stiffness initially increase (compare the n = 
0.12 and n = 0.3 dispersions in Figs. 2 and 3). Fig. 4 
shows D(n), obtained by fitting the quadratic behavior 
D{n)Q 2 to the long wavelength numerical dispersions, 
for finite exchange interaction J/t = 10. As can be 
seen in Fig. 4, the RPA predicts an initial increase of 
the spin wave stiffness with n followed by a decrease. 
However, the carrier-magnon scattering reduces the spin 
wave stiffness and changes the above concentration de- 
pendence significantly, especially in 2D and 3D (see Fig. 
4) . The difference from the RPA behavior is particularly 
striking for the 0(1/ S 2 ) contribution to D(n). The lat- 
ter is significantly suppressed as compared to the rest of 
the approximations of carrier-magnon scattering. This 
particularly strong softening indicates that the 0(1/ S 2 ) 
result significantly underestimates T c and the stability of 
the fully polarized ferromagnetic state. We note that, 
due to the variational nature of the full three-body cal- 
culation, the exact stiffness will be smaller (softer) than 
the corresponding results of Fig. 4 (solid curve). 

As the electron concentration increases, we see from 
Figs. 2 and 3 that the different approximations start to 
deviate substantially from each other. This is clear for 
the 2D system (Fig. 3), while in ID the differences de- 
velop for higher electron concentrations. Compared to 
the full three-body variational calculation ($ ^ 0, ^ 7^ 
0), the ladder and 0(1/ S 2 ) (non-variational) approxima- 
tions give softer spin wave energies, while the variational 
$ = 0,* 7^ and RPA (<f> = * = 0) wavefunctions 
give higher spin wave energies. The large differences be- 
tween the above dispersions point out the importat role 




FIG. 4: Spin wave stiffness for J = lOi as function of electron 
concentration and system dimensionality. 



of carrier-magnon correlations for such electron concen- 
trations. The difference between the full three-body cal- 
culation (or the $ = 0, 'F ^ calculation) and the ladder 
approximation in Figs. 3(c) and 3(d) shows that the 
three-body correlations are signifcant, while the differ- 
ences from the 0(1/ S 2 ) result show the importance of 
the multiple carrier-magnon scattering processes (vertex 
corrections). One can conclude from the above compar- 
isons that the Fermi sea hole-magnon scattering cannot 
be neglected for such concentrations. Furthermore, as 
can be seen in Fig. 3, the different aproximations bound 
the full three-body result. This is particularly useful for 
the 3D system, where the full three-body variational cal- 
culation could only be performed for relatively small N 
lattices with N ~ 10 (due to the large number of vari- 
ational parameters tj^). On the other hand, Fig. 5 
shows the spin wave dispersions for a rather large (50 3 ) 
3D lattice, obtained by using the RPA, 0(1/ S 2 ), and 
$ = 0, ^ 7^ approximations. By comparing the spin 
wave dispersions in Figs. 5 and 3, obtained for the same 
parameters, we see that the trends as function of n are 
qualitatively similar in the 2D and 3D systems. 

The differences between the approximations studied 
here (which are due to the correlations) are very pro- 
nounced in the electron concentration range 0.5 < n < 
0.8 for intermediate exchange interaction values, i.e. in 
the parameter regime relevant to the manganites. This 
can be seen more clearly in Fig. 6, which shows the 2D 
spin wave dispersions obtained with the different approx- 
imations along the main directions in the Brillouin zone 
for the typical parameters n = 0.7, J = 8t. Fig. 6 com- 
pares the dispersions along the directions (0, 0) — ► (tt, 0) 
(T- X), (tt, 0) -> (tt, tt) (X - M ), and along the diago- 
nal (0,0) — > (w, tt) (r — M). The discrepancies between 
the different approximations are very large along T — X 
but much smaller along the other directions. For exam- 
ple, the RPA fails completely along the direction T — X, 
where the full three-body calculation shows a striking 
spin wave softening that is most pronounced close to the 
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FIG. 5: Spin wave dispersions in 3D for J/t = 10 and electron 
concentrations similar to Figs. 2 and 3. 



X point. Such a strong effect, much stronger than the 
softening at small n, only occurs for intermediate elec- 
tron concentrations 0.4 < n < 0.7 and is underestimated 
by the <& = 0, $ ^ variational calculation. On the 
other hand, the 0(1/S 2 ) dispersion for the parameters 
of Fig. 5 shows instability at long wavelengths (negative 
stiffness) rather than softening at the zone boundary. 

To see the origin of the above spin wave softening, wc 
show in Fig. 7 the spin wave dispersion for a slightly 
smaller J/t than in Fig 6. The spin wave energy now 
becomes negative in the vicinity of the X point, while 
the magnon stiffness remains positive. This variational 
result allows us to conclude instability of the fully po- 
larized ferromagnetic state due to the X point magnons. 
The strong zone boundary softening is a precursor to this 
instability. Fig. 7 also compares the full three-body and 
RPA calculations for two different values of N with fixed 
n. For n = 0.7, our results have converged reasonably 
well even for TV ~ 10 and thus reflect the behavior in the 
thermodynamic limit. 

The above zone boundary instability occurs in 2D and 
3D for intermediate electron concentrations (0.4 < n < 
0.7 for the 2D three-body calculation), where a strong 
magnon softening and short magnon lifetimes were ob- 
served in the manganitesi 24 i 25 i 26 i 27 i 28 i 29 i 30 Although the 
other approximations can also give an instability, this 
occurs within a more limited range of n and J than for 
the full three-body calculation (discussed further below). 
The magnitude and concentration dependence of the soft- 
ening also depends on the local Hubbard repulsion (Hu) 
and direct superexchange (H super ) interactions and the 
bandstructure (to be considered elsewhere). We note 
that spin wave softening at the zone boundary of elec- 
tronic origin was obtained before within the one-orbital 
Hamiltonian for finite values of J/t by including these ad- 
ditional effectsAS The main difference here is that our 
calculation is variational (and thus allows us to draw def- 




FIG. 6: Spin wave dispersion along the different directions in 
Brillouin zone for n — 0.7, J = 8t: Comparison of the different 
approximations. 



inite conclusions by guaranteeing that the exact magnon 
energies are even softer than the calculated values) and 
our effect was obtained by using the simplest possible 
Hamiltonian (Hjj = H super = 0). 

It is important to note here that the strong magnon 
softening and instability only occur for finite values of 
J/t and disappear in the strong coupling limit J — > oo. 
This can be seen in Fig. 8, which compares the 2D 
magnon dispersions for n = 0.7 and different values of 
J/t > 1 with the result obtained by expanding Eas. JA8|) 
and (J5J in the limit J — * oo. As can be seen in Fig. 8, 
the zone boundary magnon softening disappears with in- 
creasing J/t. The magnon dispersions converge slowly 
to the strong coupling result, which is reached only for 
J/t ~ 1000. Since the typical exchange interaction val- 
ues in the manganites are of the order of J/t < 10, we 
conclude that the manganites are far from the J — > oo 
limit. Noting in Fig. 8 that the zone boundary magnon 
softening has disappeared completely for J/t > 20, we 
see that the finite J/t effects play an important role in 
the manganite spin dynamics. 

We now turn to the 3D system, where the full three- 
body variational calculation faces computational difficul- 
ties due to the large number of variational parameters 
\t5„. As can be seen in Fig. 7, in the 2D system the 
magnon dispersion results for n — 0.7 have already con- 
verged reasonably well for N ~ 10. We therefore expect 
that, in 3D, a calculation for a N x N x N lattice with 
N ~ 10 should give reasonable results. Fig. 9 show the 
3D magnon dispersions obtained this way for N = 8, 
n ~ 0.7, and J/t — 14 using the different approxima- 
tions. Fig. 9 shows similar 3D magnon behavior as in 
the 2D system (Fig. 5): magnon softening close to the 
X point and significant deviations between the differ- 
ent approximations along T — X even for this relatively 
large J/t = 14. Similar to the 2D and ID systems, the 
0(1/ S 2 ) dispersion and the carrier-localized spin scat- 
tering (<f> = 0, 7^ 0) variational results bound the full 




FIG. 7: Spin wave dispersion for a N x N 2D lattice along 
different directions (n = 0.7, J = 7.5t). Convergence with 
system size N is very good for this n. 



FIG. 8: Spin wave dispersion in 2D, obtained from the full 
three-body calculations, for n — 0.7 and increasing values of 
J/t. Convergence to the strong coupling limit is slow. 



three-body magnon dispersion. 

Next we turn to the stability of the fully polarized 
half metallic ferromagnetic: state. Our results show two 
different instabilities due to the exchange interaction 
(Hjj = H super = 0). The first is the X-point zone 
boundary instability discussed above, which occurs in 
2D and 3D for intermediate electron concentrations. For 
low or high electron concentrations (for all concentra- 
tions in ID), we find a second instability with respect 
to long wavelength spin waves (negative magnon stiff- 
ness) similar to previous calculations. In this case, the 
minimum magnon energy occurs at a finite momentum 
value instead of Q = 0. This momentum increases with 
n and becomes ir at n = 1 (antiferromagnetic order at 
half filling). This result implies instability to a spiral 
state, while the system can further lower its energy by 
phase separatingi^S Due to the variational nature of our 
calculation, it is guaranteed that, if the magnon energy 
becomes negative for J = J c (n), the ground state of the 
Hamiltonian H for all J < J c (n) is not the half metallic 
state \F). The phase diagrams of Fig. 10 describe the 
stability of this state against spin wave excitations. The 
most striking feature in Fig. 10 is the large shift (in- 
crease) of the ferromagnetic phase boundary, J c (n), as 
compared to the RPA, due to the carrier-magnon scat- 
tering. Furthermore, the different approximations of the 
carrier-magnon scattering lead to significant differences 
in J c {n). By comparing the shape of J c (n) between the 
ID and 2D/3D systems, we see that, in the latter case, 
J c (n) develops a plateau-like shape within an intermedi- 
ate concentration regime (see Fig. 10(d)). This feature is 
absent in the ID system, where there is no zone bound- 
ary instability. This plateau occurs for 0.4 < n < 0.7 in 
2D (full three-body calculation) and for 0.25 < n < 0.6 
in 3D (<£> = 0, ^> ^ three-body calculation). It is much 
less pronounced for the 0(1/ S 2 ) and RPA calculations 

For small n, J c (n) is small, implying enhanced stabil- 



ity of the ferromagnetic state in the concentration regime 
relevant, e.g., to III-Mn-V semiconductors!^ This stabil- 
ity is a remnant of the fact that, in the exactly solvabe 
limit N e = 1, the ferromagnetic state \F) is the ground 
state for all values of J/t. J c (n) increases more slowly 
with n in ID than in 2D and 3D. This implies enhanced 
stability of the fully polarized ferromagnetic state, partly 
due to the lack of zone boundary instability in ID. Fig. 
10(d) shows the 2D phase diagrams for the electron con- 
centrations relevant to the manganites. The full three- 
body variational calculation gives J c (n) ~ 7 — 8t in this 
regime, close to the high end of the values quoted in the 
literature^ Therefore, the simple double exchange Hamil- 
tonian predicts that the manganites lie in a regime that is 
close to the instability of the ferromagnetic state. In this 
regime, the correlations, vertex corrections, and finite J 
effects play an important role in the spin dynamics. 

Fig. 10 also compares the phase diagrams obtained 
by using the different approximations implemented here. 
The 0(1/ S 2 ) calculation underestimates the stability of 
the fully polarized ferromagnetic state, while the RPA 
and carrier-localized spin scattering ($ = 0, 'J ^ 0) vari- 
ational results overstimate the stability. For low con- 
centrations n < 0.3, all the different treatments of the 
carrier-magnon scattering predict a similar ferromag- 
netic phase bounday. It is clear from Fig. 10(d) that, in 
the electron concentration regime 0.5 < n < 0.8 relevant 
to the manganites, the RPA significantly overestimates 
the stability of the ferromagnetic state. For example, 
for n ~ 0.5, the RPA underestimates J c (n) by 100% as 
compared to the full three-body variational calculation. 
Finally, close to half filling n = 1, the two variational 
results give magnon energies similar to the RPA, which 
becomes exact for n = 1. On the other hand, the 0(1/ S 2 ) 
approximation fails in this high concentration regime. 

Finally we discuss the relevance of our calculation to 
the spin wave dispersion observed experimentally in the 
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FIG. 9: Spin wave dispersion for a N x N x N 3D lattice 
along different directions (n = 0.7, J = 14f, and N = 8). 
r = (0, 0, 0), X = (tt, 0, 0), M = (tt, tt, 0), J? = (tt, tt, tt). 

quasi-2D and 3D manganites. The experimental re- 
sults are typically analyzed by fitting the short range 
Heisenberg dispersion to the long wavelength experi- 
mental dispersion and then comparing the two close to 
the zone boundary, 23 i 24 i 25 . 26 i 27 i 28 i 29 i 30 i 38 This compari- 
son showed that the Heisenberg model fails to describe 
the experimental results in the overdoped manganites 
(typically for 0.5 < n < 0.7), but fits well in the un- 
derdoped samples (n > 0.7). This failure is due to 
the strong magnon softening close to the zone bound- 
ary (X-point)j 24 i 25 i 26 i 27 i 28 i 29 ' 30 whose physical origin is 
currently under debatei 24 i 25 i 26 i 27 i 28 i 29 i 30 i 31 i 32 i 33 Here we 
compare our numerical results with the Heisenberg dis- 
persion Wq cis , obtained by fitting to the long wave- 
length numerical results, by introducing the parameter 
A = bj^ els /u)x — 1) where (jj^ els and u>x are the Heisen- 
berg and numerical magnon energies calculated at the 
X-point. |A| thus measures the magnitude of the devi- 
ations from Heisenberg behavior at the zone boundary. 
For example, |A| ~ 1 means 100% deviation, A > 
means magnon softening at the zone boundary, as com- 
pared to the Heisenberg dispersion with the same stiff- 
ness, while A < implies zone boundary hardening. Fig. 
11 compares A(n) obtained from our different approxi- 
mations. With the exception of small values of J/t, the 
RPA gives small deviations from Heisenberg behavior, 
mostly a hardening at the zone boundary ( A < 0, see 
Fig. 10), and predicts a weak concentration dependence 
of A(n). This similarity between the RPA and Heisen- 
berg dispersions is expected for J 3> t since the two co- 
incide in the strong coupling limit J — > oo (see Eq.@). 

The magnon-electron scattering leads to larger devi- 
ations from Heisenberg ferromagnet spin dynamics and 
enhances A(n) (see Figs. 11(a) and 11(b) for 2D and 3D 
respectively). In order to compare with the experiment, 
the value of J/t must be chosen so that \F) is stable up 
to n ~ 0.8, where a metallic ferromagnetic state is ob- 



FIG. 10: Phase diagram due to the spin wave instability and 
comparison between the different approximations discussed in 
the text, (a) ID system, (b) 2D system, (c) 3D system, and 
(d) 2D system in the electron concentration range relevant to 
the manganites. 



served experimentally. For J/t ~ 10, this is the case for 
the full three body calculation, while larger values of J/t 
are required to achieve stability for n ~ 0.8 with respect 
to the 0(1/ S 2 ) magnons. Figs. 11(a) and 11(b) compare 
the behavior of A(n) for the different approximations in 
the 2D and 3D systems respectively. The 0(1/ S 2 ) calcu- 
lation gives magnon hardening rather than softening in 
the concentration range of interest, similar to the strong 
coupling results of Ref. |l(J. This is in contrast to A(n) 
obtained by using the full three-body calculation, shown 
in Fig. 11(a) for the 2D system. In this case, the magnon 
hardening for n < 0.5 (A < 0) changes to magnon soft- 
ening for 0.5 < n < 0.7 (A > 0) and then back to a 
small magnon hardening for n > 0.7. This behavior with 
n is consistent with the experimental trends. Although 
magnon softening at the X point can be obtained us- 
ing other approximations, the full three-body calculation 
gives such an enhanced effect within the range of inter- 
mediate electron concentrations of interest and for values 
of J/t such that the fully polarized ferromagnetic state is 
stable for 0.5 < n < 0.8 (where it is observed experimen- 
tally). The above behavior of A(n) is not reproduced in 
the strong coupling limit J — > oo, where magnon hard- 
ening is obtained. It arises from the interplay of the 
X-point instability and the plateau-like shape of J c (n), 
Fig. 10, induced by the correlations. On the other hand, 
for J — lOi, the carrier-localized spin scattering approx- 
imation ($ = 0, f / variational wavefunction) gives 
A(n) that, more or less, follows the RPA behavior (see 
Figs. 11(a) and 11(b)). As J/t decreases, magnon soften- 
ing, A > 0, can also be obtained with this approximation 
over a range of electron concentrations in both 2D and 
3D (see Figs. 11(c) (2D) and 11(d) (3D)). However, for 
such J/t, the ferromagnetic state is unstable for n > 0.6, 
i.e. in a regime where ferromagnetism is observed exper- 
imentally. We expect that the precise behavior of A(n) 



11 






1 

(b) 3D, j =10t 








\ "V 
\ \ 


3-body 


■ ■ ■ RPA 


\ 


— - 4>=0 




■ - 0(1/S 2 ) 


\ 

1 



0.1 

0.3 



1 1 


' (d) 3D, J =6t 




- — — \ 




\ 


/ 


\ 


/ 


\ ~ 


./ 


' • . ,\ 


i 


V 

1 , 1 



6 -<H 



FIG. 11: Deviation A(n), defined in the text, from Heisenberg 
ferromagnet spin dynamics in the 2D and 3D systems for fixed 
J/t — 10 and different electron concentrations. 



in the realistic materials will also depend on Hjj , H super , 
and the bandstructure effects (to be studied elsewhere). 
Here we point out that at least three-body correlations 
must be included for a meaningful comparison to the ex- 
periment. 



IV. CONCLUSIONS 

In this paper we presented a nonperturbative varia- 
tional calculation of the effects of magnon-Fermi sea pair 
correlations on the spin wave dispersion for the simplest 
possible double exchange Hamiltonian. Our theory treats 
exactly all three-body long range correlations between 
an electron, a Fermi sea hole, and a magnon excitation. 
We achieved this by using the most general variational 
wavefunction that includes up to one Fermi sea pair exci- 
tations. Since the contribution of multipair Fermi sea ex- 
citations is suppressed by powers of 1/S, one could alter- 
natively think of our calculation as putting the 0(1/ S 2 ) 
result, which treats the one Fermi sea pair contribution 
perturbatively within the Born approximation, on a vari- 
ational nonperturbative basis. Our theory (i) becomes 
exact in the two limits of one and N d electrons and should 
therefore interpolate well between the low concentration 
and half filling limits, (ii) converges well with system 
size and thus applies to the thermodynamic limit, (iii) 
becomes exact in the atomic limit (i=0), conserves mo- 
mentum exactly, and treats both short and long range 
correlations on equal basis; it should therefore interpo- 
late well between the strong and weak coupling limits, 
which is important given the relatively small values of 
J/t in the manganites, and (iv) contains the well known 
0(1/ S 2 ) and RPA results as limiting cases. In this pa- 
per we studied, among others, (i) the shape of the spin 
wave dispersion and ferromagnetic phase boundary for 
different system dimensionalities (ID, 2D, and 3D), (ii) 
the deviations from the strong coupling double exchange 



limit, and (iii) the role of up to three-body correlations 
and nonperturbative vertex corrections on the spin dy- 
namics. By comparing the full three-body variational 
calculation to a number of approximations (RPA, 1/S 
expansion, ladder diagram treatment of two-body cor- 
relations, and carrier-localized spin rather than carrier- 
magnon scattering), we showed that the correlations play 
an important role on the spin excitation spectrum, the 
stability of the ferromagnetism, and the shape of the fer- 
romagnetic phase boundary in the parameter range rel- 
evant to the manganites. Importantly, the correlations 
lead to spin dynamics that differs strongly from that of 
the short range Heisenberg ferromagnet for intermediate 
electron concentrations. 

Our main results can be summarized as follows. First, 
the different approximations lead to substantial differ- 
ences in the spin wave dispersion and ferromagnetic phase 
boundary for electron concentrations above n ~ 0.3 and 
intermediate values of J/t, which includes the parameter 
range relevant to the manganites. These large differences 
come from the correlations, which cannot be neglected, 
and imply that variational calculations should be used 
if possible in order to draw definite conclusions. Sec- 
ond, we find that, depending on n, there are two possi- 
bile instabilities of the ferromagnetic state toward spin 
wave excitations: instability driven by a negative spin 
stiffness and instability at large momenta, close to the 
A-point zone boundary, with positive stiffness. The lat- 
ter instability only occurs in the 2D and 3D systems, 
for electron concentrations n < 0.7 and finite values of 
J. The three-body carrier-magnon correlations enhance 
this effect. As a precursor to the above zone boundary 
instabilty, we find a strong magnon softening at the X- 
point, which should be accompanied by short magnon 
lifetimes. Third, by comparing to the Heisenberg disper- 
sion obtained by fitting to the long wavelength numerical 
results, we find strong deviations from the spin dynam- 
ics of the short range Heisenberg model. By choosing the 
exchange interaction so that the fully polarized ferromag- 
netic state is stable up to n > 0.8 as in the experiment, we 
show that the full three-body 2D calculation gives strong 
magnon softening at the A point for 0.5 < n < 0.7, which 
changes into a small hardening for n > 0.7. This is sim- 
ilar to the behavior observed in the manganites. Our 
work provides new insight into the spin dynamics in the 
manganites and can be extended to treat related ferro- 
magnetic systems (such as e.g. the III-Mn-V magnetic 
semiconductors) that are far from the double exchange 
strong coupling limit. Our calculations imply that the 
metallic ferromagnetic state in the manganites should 
be viewed as a strongly correlated state. Finally, the 
carrier-magnon correlations studied here can also play 
an important role in the ultrafast relaxation dynamics of 
itinerant ferromagnetic systems, which is beginning to be 
explored by using ultrafast magneto-optical pump-probe 
spectroscopy4i*22i2ii 

We thank N. Papanicolaou for very useful discussions. 
This work was supported by the EU Research Training 



12 



Network HYTEC (HPRN-CT-2002-00315). 



APPENDIX A 

In this Appendix we present the variational equations 
that determine the wavefunction amplitudes X", W^, 
and <J>Q. These are obtained by minimizing the vari- 
ational energy (F\MqHMq\F), where Mq is given by 
Eq.Q, with respect to the above variational amplitudes. 
The normalization condition (F\MqMq\F) — 1 is en- 
forced via a Lagrange multiplier, which coincides with 
the variational magnon energy wq. Similar to the three- 
body Fadeev equations, the resulting variational equa- 
tions are equivalent to solving the Schrodinger equation 
within the subspace spanned by the states SqI-F) and 

CQ + „jC„-f \F), which describe all possible configurations 
with one reversed spin and no Fermi sea excitations, and 
the magnon-Fermi sea pair states c' a ^c u -\ SQ +l/ _ a \F) and 

clT^TCQ+M+tz-ai^Tl-F) to which a magnon or spin flip 
excitation can scatter with the simultaneous excitation 
of a Fermi sea pair. We note that the above momentum 
space basis ensures the conservation of momentum and 
total spin. The explicit form of the variational equations 
is obtained after straighforward algebra by projecting the 
Schrodinger equation [H, Mq]\F) = ojq\F) in the above 
basis after calculating the commutator [H , Mq] by us- 
ing Eq.Q for Mq and noting that H\F) = ( we take 
the energy of \F) as zero). The variational equation that 
gives the energy ujq reads 



Jn 

T 



j 

2N 



Yx^ + — Y^ Q 



(Al) 



The last term in the above equation describes the contri- 
bution due to the carrier-magnon scattering. The first 
two terms on the rhs give the RPA magnon energy if X® 
is substituted by its RPA value, obtained for = 0. 
The carrier-magnon scattering renormalizes X® as com- 
pared to the RPA result: 



(JS + e„ +Q -e v - w q ) X® = JS 



(A2) 



The first term on the rhs of the above equation gives 
the RPA contribution to X™, while the second term, as 
well as the correlation contribution to the magnon energy 
luq , describe the effects of the magnon-carrier scattering. 
The RPA is obtained from Eqs. I|A1|) and (|A2|> after 
setting * Q = $ Q = 0: 

< A = ^-^Y, X ? RPA > ( A3 ) 



XQBPA = 



2N ■ 
JS 



JS + c„+Q -£ u - Uq PA ' 



(A4) 



This full RPA result can also be obtained as the zeroth or- 
der contribution to an expansion in powers of 1/L, where 
L is the number of electron flavors and corresponding 
degenerate electron bands42i 

The scattering amplitude x f ,( ^ is determined by the 
variational equation 



Jn 



e„-e a )*% = ^(l-X?) 



--yf« - -Ly^,- J-Yqq, . ( A 5) 



The first term on the rhs of the above equation gives 
the Born scattering approximation contribution to the 
carrier-magnon scattering amplitude, which is the only 
one that contributes to 0(1/S 2 ). The next two terms 
describe the effects of the multiple electron-magnon (sec- 
ond term) and hole magnon (third term) scattering. Fi- 
nally, the last term comes from the electronic contribu- 
tion to the scattered magnon, i.e. from the coherent ex- 
citation of a spin-| electron to the spin-| band. The am- 
plitude $^ of the latter contribution to Eq.Q) is given 
by the variational equation 

(JS + EQ+^+^-a - e M + e a — e v — ujq) = 

JS{9* (A6) 

We note that, in the strong coupling limit J — > oo, 

^aufj. ~ and tlie last tw0 terms in E q- © de- 

scribe the scattering of a Fermi sea pair with the strong 

coupling RPA magnon created by the total spin lower- 
ing operator Sq + -7= J2v c Q+vi c »1- Our general wave- 
function Eq. J3J does not assume an RPA magnon and 
includes corrections to the strong coupling limit that are 
important for the values of J/t relevant to the mangan- 
ites. 

A closed equation for the carrier-magnon scattering 
amplitude can be obtained by substituting in Eq. 
151) the expressions for <j>Q and X® obtained from Eqs. 
03) and (|A2|) . Defining the excitation energy 



E 



EQ+v'+v-c 



'av — W Q + £ v — 



J 

2iV *—f JS + £q+v'+v-u — £v' + s a — e v — wq 
we thus obtain the following equation: 



(A7) 



J e v+Q — £v — Wq 



av av 



2N JS + £j,+q - £ v — wq 



2N ^ av 



£ Q+iH 



v' — a &v' +e a — e 



a' 



JS + £Q+v+v'-a — £v> + £ a — S u — LOq 



(A8) 



The above equation describes up to three-body correla- 
tions between the magnon and a Fermi sea pair. The 
first term on the rhs describes the bare carrier-magnon 
scattering amplitude. This is renormalized by the mul- 
tiple scattering of a Fermi sea electron (second term on 
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the rhs) and a Fermi sea hole (last term on the rhs). 
These two contributions describe vertex corrections of 
the carrier-magnon interaction. Eas. (|Al|) and l)A8|l were 
solved iteratively until convergence for the spin wave en- 
ergy was reached. 



APPENDIX B 

In this Appendix we identify the three-body correla- 
tion contribution to the carrier-magnon scattering ampli- 
tude W***, Eq. (| AS|) - and distinguish it from the two-body 
multiple scattering contributions. 

We note from Eq. 1A8|I . that \E' Q has the form 



The first factor on the rhs of the above equation comes 
from the electron-magnon two-body ladder diagrams 
summed to infinite order. The last term in the second 
factor describes the coupling of the electron-magnon and 
hole-magnon scattering channels. This coupling comes 
from the three-body correlations. An analogous equa- 
tion for T h can be obtained by substituting Eq. (|B1|) into 
Eq. I|B2|) . In the case of the simpler variational wavefunc- 
tion $ = 0, ^ 7^ 0, which describes the carrier-localized 
spin scattering contribution, the calculation of T h simpli- 
fies by noting from Eq. l|A5l) and the definition Eq. ljBll) 
that = r^. The corresponding variational equation 
can be obtained by setting = in Eq. (|A5I) : 



J 



pe pn 

u av 



av 2NA%JS + e v+Q -e v 



A Q ' 
(Bl) 

where we introduced the electron vertex function Eq.© 
and the hole vertex function 



av 9 AT 2-^1 av 



2N 



JS + £Q+v+v' -a — £V + £q — S v — UlQ 

(B2) 

Substituting Eq. (jBTJ into Eq. © we obtain after some 
algebra that 



V 



_J_ 

2N 



E 



AN 2 *—f JS + £„'+q — e v i — ujq A*^ , 2N^ A Q 



2N JS + e u+Q - e„ - w Q A Q , 



E A Q 



AN 2 \JS + e v+Q -e v -uci) A^ 



J 



a' ol'v 



2N JS + £„ +Q - e u — ojq 2—f A 



(B3) 



The first factor on the rhs comes from the hole-magnon 
ladder diagrams, while the coupling to T e comes from the 
three body correlations. 
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